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ABSTRACT 

We present metallicity distribution functions (MDFs) for the central regions of eight dwarf satellite 
galaxies of the Milky Way: Fornax, Leo I and II, Sculptor, Sextans, Draco, Canes Venatici I, and 
Ursa Minor. We use the published catalog of abundance measurements from the previous paper in 
this series. The measurements are based on spectral synthesis of iron absorption lines. For each 
MDF, we determine maximum likelihood fits for Leaky Box, Pre-Enriched, and Extra Gas (wherein 
the gas supply available for star formation increases before it decreases to zero) analytic models of 
chemical evolution. Although the models are too simplistic to describe any MDF in detail, a Leaky 
Box starting from zero metallicity gas fits none of the galaxies except Canes Venatici I well. The 
MDFs of some galaxies, particularly the more luminous ones, strongly prefer the Extra Gas Model 
to the other models. Only for Canes Venatici I does the Pre-Enriched Model fit significantly better 
than the Extra Gas Model. The best-fit effective yields of the less luminous half of our galaxy sample 
do not exceed 0.02 Z®, indicating that gas outflow is important in the chemical evolution of the less 
luminous galaxies. We surmise that the ratio of the importance of gas infall to gas outflow increases 
with galaxy luminosity. Strong correlations of average [Fe/H] and metallicity spread with luminosity 
support this hypothesis. 

Subject headings: galaxies: dwarf — galaxies: abundances — galaxies: evolution — Local Group 



1. INTRODUCTION 

The star formation history of a galaxy shapes the 
metallicity distribution of its stars. Therefore, simply 
counting the number of stars in different bins of metal- 
licity in a galaxy is a way to quantify the gas dynam- 
ics during the history of star formation in the galaxy. 
How much gas was accreted by gravitational attraction? 
How much gas left the galaxy from supern ova winds (e.g., 
iDekel fc Silk 1986: lGovernato et al.ll2010l) or tidal or ram 
press ure stripping (|Lin Sz Faberl 119831 : IMarcolin i et al.l 
2003) from interaction with the Milky Way (MW)? 

The most basic approach to answering these questions 
is to fit an analytic model of chemical evolution to the ob- 
served metallicity distribution function (MDF). For ex- 
ampl e, one could assume th at the galaxy is a "closed 
box" ([Talbot fc Arnett]|197lT) . In other words, the galaxy 
begins its life with a fixed amount of gas. It loses gas 
only to t he formation of stars, and it does not acq uire 
new gas. IVan den Berghl (IT962I) and[Sclrmidj (fl963f) fa- 
mously applied this model to the metallicity distribution 
of G dwarfs in the solar neighborhood to find that the 
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MW disk is not a closed box. Instead, it e xperiences 
more complicated gas dynamics. iPagell (j!997t ) described 
some more complex analytic models, including some that 
incorporate the accretion of external gas during the life- 
time of star formation. 

Our aim is to examine the metallicity distributions of 
dwarf spheroidal satellite galaxies (dSphs) of the MW 
in order to reveal how gas infall and outflow affected 
their star formation histories. We use the [Fe/H] mea- 
surement s for individual star s in eight MW dSphs from 
iPaper III (jKirbv et al.1 l2010b[ ) . The measurements are 
based on spectral synt hesis of iron lines from medium 
resolut ion spectra. IKirbv. Guhathakurta. &: Sneder] 
( 2008a) described the measurements in detail. We mod- 
ified th e procedure for determin ing [Fe/H] in Papers I 
and II (IKirbv et al.ll2f)09l I2010H) . 

Several previous spectroscopic studies have examined 
MW dSph metallicity distributions in the context of 
star formation history. Most of them rely on a calibra- 
tion between the summed equiva lent widths of the C a 
infrared triplet and [Fe/H] (e.g., iRutledge et al.lll997ft . 
iTolstov et al.l (2001) first published a significant sample 
of Ca triplet metallicities for the Sculptor and Fornax 
dSphs. From the widths of the metallicity distributions 
and the age spreads apparent from their broadband col- 
ors combined with spectroscopic metallicities, they con- 
cluded that both dSphs experienced extended star for- 
mation, unlike globular clusters. Since then, the Dwarf 
Abundances and Radial Velocities Team (DART) have 
measured MDFs for Carina and Sextans in addition to 
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vations with a multi-object spectrograph. They found 
extre mely metal-poor stars i n the ultra-faint dSph sam- 
ple of lSimon fc Gehal (120071) . IShetrone et all (12001) also 
applied spectral synthesis to medium-resolution spectra. 
They obt ained [Ee / H] measurements for 27 red giants in 
Leo II. In iPaper II (jKirbv et al.1 120091 ) . we presented the 
MDF for the Sculptor dwarf galaxy based on the spectral 
synthesis of iron lines for 388 red giants. 

In this paper, we extend our analysis of Sculptor to 
seven additional MW dSphs: Fornax, Leo I and II, Sex- 
tans, Draco, Canes Venatici I, and Ursa Minor. The 
sample of 2961 stars in eight galaxies permits a com- 
parative look at star formation histories, particularly as 
the properties of the MDFs change with dSph luminos- 
ity. Furthermore, we examine how the distributions of 
[Fe/H] change as a function of distance from the center 
of each dSph. 

We describe three analytic chemical evolution models 
in Sec. [H They are Pristine, Pre-Enriched, and Extra 
Gas Models. In Sec. [3l we analyze the MDF of each 
dSph in detail. We discuss how each chemical evolu- 
tion model may or may not apply to different dSphs. 
The results of previously pub l ished numerical models 
(jLanfranchi &: Matteuccil I2003L 120041 ) are compared to 
the observed MDFs in Sec. |H Section [S] is devoted to 
radial gradients and their relevance to the dwarf galax- 
ies' star formation histories. We quantify trends of the 
MDF properties with luminosity in Sec. |6] and the trends 
of the most likely chemical evolution model parameters 
with dSph properties in Sec. [7] Finally, in Sec. [H we 
summarize our conclusions, point out shortcomings in 
the chemical evolution models and our conclusions from 
them, and suggest how our work may be improved in the 
future. 

2. ANALYTIC CHEMICAL EVOLUTION MODELS 

The theory of galactic chemical evolu tion has pro- 
gresse d significantly since iTinslevI (Il980l) codified the 
field. lLanfranchi fe Matteuccil ( 20041) applied sophisti- 
cated multi-element models to the then scant spectro- 
scopic stel lar abundance measurements in MW satellite 
galaxies. iMarcolini et al.1 (|2006l 120081 ) created three- 
dimensional, hydrodynamic models of isolated dwarf 
galaxies, but the sample size of observatio ns was inad- 
equate to test their models. Only recently, iRevaz et al.l 
(2009) modeled the abundance distributions of several 
MW satellites. They obtained a reasonable agreement 
between their predictions and the large samples of the 
Dwarf Abundances and Radial Velocities Team's abun- 
danc e measur ements in Fornax and Sculptor. 

In IPaper II of this series, we fit two different chemi- 
cal evolution models to the MDF of the Sculptor dwarf 
galaxy: a Simple Model and an Extra Gas Model. The 
Simple Model represented a leaky box. Although the 
gas was allowed to leave the galaxy, the galaxy never ac- 
quired new gas. Gas outflows do not affect the functional 
form of the MDF as long as the nucleosynthetic yield p is 
assumed to be the effective yield and not the true metal 
yield of the stars. In the Simple Model, the initial gas of 
the galaxy was also allowed to be be pre-enriched with a 
metallicity of [Fe/H] . 

We distinguish between the Simple Model with pristine 
initial gas ([Fe/Hl n = ~oo) an d pre-enriched initial gas 
([Fe/H] is finite). iPagel (fl997h gave the functional form 



of the MDF of the Pre-Enriched Model: 
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where p is in units of the solar metal fraction (Zq). For 
the Pristine Model, one term vanishes: 
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In both the Pristine and Pre-Enriched Models, the 
peak of the MDF increases with p. However, a metal-rich 
dwarf galaxy did not necessarily host supernovae with 
higher yields than a more metal-poor galaxy, even if the 
Pristine Model is a good description for both galaxies. 
Because p represents the effective yield, it encapsulates 
both the supernova yield and gas outflow. Larger su- 
pernova yield increases p, and more intense gas outflow 
decreases p. Therefore, an equally valid interpretation of 
the hypothetical MDFs is that the metal-rich (large p) 
galaxy retained gas more effectively than the metal-poor 
galaxy. 

The Extra Gas M odel is the Best A ccretion Model of 
iLvnden-Bell (| 19751 also see lPagellll997l ). Unlike the Pris- 
tine and Pre-Enriched Models, it allows the galaxy to 
access an additional supply of gas available for forming 
stars during or between other episodes of star formation. 
For simplicity, we assume that the gas is metal-free, and 
we assume that the rate of gas infusion decays over time. 
Furthermore, we assume a relation between the gas mass 
and the stellar mass that permits an anal ytic solution 
to the differential metallicity distribution. iLvnden-Belll 
generated such a relation for which the gas mass reached 
a maximum and for which the stellar mass rose asymp- 
totically to its final value. Thes e qualitative ch aracteris- 
tics matched the simulations of iLarsonl ()1974al ). In this 
model, the gas mass g in units of the initial gas mass is 
related quadratically to the stellar mass s in the same 
units: 



(3) 



where M is a parameter greater than 1. In the special 
case where all of the gas is converted into stars and p is 
the true yield, M equals the final stellar mass in units 
of the initial mass of gas at the onset of star formation. 
When M = 1, Eq.[3] reduces to g = 1 — s, which describes 
a closed box, wherein the gas mass is depleted only by 
star formation. Therefore, the Extra Gas Model reduces 
to the Pristine Model when M = 1. Otherwise, an in- 
crease in M represents an increase in the amount of gas 
the galaxy gains. The following two equations describe 
the differential metallicity distribution: 
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Equation 2] is transcendental, and it must be solved for 
s numerically. The solution to s may then be put into 
E quation [5] 

Lvndcn-Bcll (1975) named his model the Best Accre- 
tion Model because he assumed that new gas became 
available for star formation by the accretion or infall of 
gas onto the galaxy. We have chosen the more general 
name Extra Gas Model. The model is too simplistic to 
distinguish between different mechanisms of cold gas in- 
fusion. For example, the galaxy could contain hot gas. If 
that gas cools between episodes of star formation, then 
it becomes available to form stars. However, the Extra 
Gas Model does require the newly available gas supply 
to be metal-free, or at least much more metal-poor than 
the bulk stellar metallicity at any given time. 

Allowing for extra gas complicates the interpretation of 
the peak of the MDF. Both a larger supernova yield and 
smaller increase in the supply of pristine gas would in- 
crease the peak of the MDF. However, gas outflow would 
still decrease the peak of the MDF. Thankfully, under 
the assumption of instantaneous mixing, the parameter 
M uniquely quantifies the amount of extra gas, leaving 
p to be degenerate only between the supernova yield and 
gas outflow. 

The Pristine, Pre-Enriched, and Extra Gas models all 
assume the instantaneous recycling approximation (IRA) 
and the instantaneous mixing approximation (IMA). The 
IRA poorly reproduces the distribution of secondary nu- 
clides, which are produced on longer timescales than pri- 
mary nuclides. Unfortunately, iron is a secondary nu- 
clide, but it is the most precisely measured of any element 
in stellar spectroscopy because it has a large number of 
absorption lines in the visible sp ectrum. The IMA ma 
not be appropriate for dSphs. iMarcolini et al.l (200 
showed that inhomogeneous pollution from Type la SNe 
affects the modeled MDF of a Draco-like dSph within 
two core radii. However, analytic forms of the differen- 
tial metallicity distributions require the assumption of 
both the IRA and the IMA. 

3. METALLICITY DISTRIBUTIONS 

We fit the Pristine, Pre-Enriched, and Extra Gas Mod- 
els to th e MDF o f each of the eight dSphs in the cata- 
log from iPaper iH . The model parameters — one for the 
Pristine Model and two for each of the Pre-Enriched and 
Extra Gas Models — were determined by maximum likeli- 
hood. Each analytic metallicity distribution was treated 
as a probability distribution (dP/d[Fe/H]) which is nor- 
malized as dP/d[Fe/B] d[Fc/H] = 1. The functional 
forms of the probability distributions were identical to 
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maximizing the likelihood function L. 



n 



exp 



dP 



3 d[Fe/H] V2^<5[Fe/H] 
([Fe/H]-[Fe/HU 2 \ 



2(5[Fe/R] i Y 



d[Fe/H] (6) 



The index i represents each star in the observed MDF. 
For computational simplicity, the most likely parame- 
ters were actually determined by minimizing the nega- 
tive, logarithmic likelihood function L. 
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We initially estimated the model parameters that max- 
imized likelihood using the Powell optimization method. 
Then we used a Metropolis-Hastings Markov Chain 
Monte Carlo algorithm to refine the fit and to estimate 
measurement uncertainties. The proposal distributions 
were normally distributed with a = 0.01 for the effective 
yield parameters and u = 0.1 for the [Fe/H] and M pa- 
rameters. We conducted 10 3 trials of the one-parameter 
Pristine Model after 10 2 burn-in trials. We conducted 
10 5 trials of the two-parameter Pre-Enriched and Extra 
Gas Models after 10 3 burn-in trials. The fiducial best- 
fitting parameters were taken to be the median values 
of all of the trials. Finally, we computed the two-sided 
68.3% confidence interval. The upper error bar was the 
value that included 68.3% of the trials above the median. 
The lower error bar was the value that included 68.3% 
of the trials below the median. 

The relative goodness of fit of one model over another 
may be quantified by the ratios of the maximum like- 
lihoods. Once the parameters for one model are de- 
termined by maximizing L, then the model that bet- 
ter describes the data will have a larger maximum like- 
lihood, i m ax- In the following sections, we quantify 
the relative goodness of fit between two models as the 
logarithm of the ratio of their maximum likelihoods: 
lnL max (Model A)/L max (Model B). 

The Pristine Model consists of only one free param- 
eter, the effective yield. Both the Pre-Enriched and 
Extra Gas Models depend on two free parameters. In 
fact, the Pristine Model is a special case of both mod- 
els ([Fe/H] Q = — oo for the Pre-Enriched Model and 
M = 1 for the Extra Gas Model). Therefore, the 
Pristine Model will never be more likely than the Pre- 
Enriched or the Extra Gas Models. (Sculptor is an ex- 
ception because we used the upper limit for [Fe/H] in 
the calculation of L max for the Pre-Enriched Model.) 
However, the ratio In L max (Pre-Enriched or Extra Gas)/ 
Lmax (Pristine) may be close to zero, indicating that the 
extra free parameter does not add significantly to the 
description of the MDF. 

Figure [T] shows the observed metallicity distributions 
of each dSph with Poisson error bars. The histograms 
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Fig. 1. — The differential metallicity distribution in each dSph expressed as a fraction of the total number of observed stars. The error 
bars represent Poisson counting statistics. The panels are arranged from left to right and then top to bottom in decreasing order of dSph 
luminosity. The black histograms show only stars with estimated uncertainties of <5[Fe/H] < 0.5. The number of such stars is given in the 
upper right corner of each panel. The blue, red, and green curves are the maximum likelihood fits to galactic chemical evolution models 
(Eqs.[T][2] and[5]l convolved with the measurement uncertainties. The dotted orange lines in some panels show predictions from numerical 
models (Sec. [2) convolved with the measurement uncertainties. 
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include measurements with estimated uncertainties on 
[Fe/H] less than 0.5 dex. Stars with larger uncertainties 
are excluded from the figures for clarity but not from the 
maximum likelihood fits. This restriction excludes fewer 
than 10 stars from the plots of the four most luminous 
dSphs. The lower stellar density in the four least lumi- 
nous dSphs necessitated targeting fainter stars. There- 
fore, the fraction of stars excluded from the plots is be- 
tween 9% and 12% in Sextans, Draco, Canes Venatici I, 
and Ursa Minor. 

In order to approximate the widening caused by mea- 
surement error, the most likely analytic metallicity distri- 
butions have been smoothed. The smoothing kernel was 
a sum of N Gaussians, where N is the number of stars 
that passed the uncertainty cut. The width (A [Fe/H]) 
of the i th Gaussian corresponded to the estimated uncer- 
tainty on the i th measurement of [Fe/H]. The smoothing 
kernel was normalized to preserve the area under the 
MDF. The kernel was constant with [Fe/H] because we 
have not observed a significant correlation between our 
estimates of <5[Fe/H] and [Fe/H] except at the very lowest 
metallicities ([Fe/H] < -3). 

Our conclusions are valid for the subset of the stellar 
populations observed. Most dSphs have radial metallic- 
ity gradients (see Sec. O such that the outermost stars 
are more metal-poor than the innermost stars. Our spec- 
troscopic observations were centered on the dSphs in 
order to maximize the number of member stars. Con- 
sequently, our samples are more metal-rich and proba- 
bly younger than the dSphs' entire stellar populations. 
This effect is especially pertinent because some dSphs 
have been known to lose their outermost , oldest stars via 
tidal s tripping by the MW (e.g., Carina, Maicwsk i et al.l 
2000tJ). Consequently, our results are not applicable to 
the en tire star formatio n histories of the dSphs, such as 
Leo I ( So hn et al.ll2007t ). that may have shed significant 
fractions of their older, more metal-poor stars. 

3.1. Fornax 

Fornax is the most luminous of the dSphs that we con- 
sider, and it the most luminous intact dSph that orbits 
the MW. Sagittarius is more luminous, but the MW 
has tidally disrupted it (llbata et al.l 119941), and it may 
have been too luminous (iNiederste-Ostholt et al.l 1 20101 ) 
and too disky (jPenarrubia et al.ll2010f ) to belong to the 
same class of galaxies as the surviving dwarf spheroidals. 
In agreement with the metallicity- luminosity relation fo r 
MW dwarf satellites ()Mateol fl998t iKirbv et al.ll2008bD . 
Fornax also shows the highest peak [Fe/H] of all of the 
MDFs in Fig.Q] It has a median [Fe/H] of -1.01. 

The Pristine and Pre-Enriched Models of chemical 
evolution poorly match the observed MDF. They do 
not allow a narrow peak, nor do they allow a sharp 
change in slope on the metal-poor end of the peak. In- 
stead, the Extra Gas Model matches the observed MDF 
much better. The logarithm of the ratios of the maxi- 
mum likelihood for the Extra Gas Model to the maxi- 
mum likelihoods for the Pristine and Pre-Enriched mod- 
els are lnI/ max (Extra Gas) /L max (Pristine) = 125.62 and 
In L max (Extra Gas) /L max (Pre-Enriched) = 75.64. The 
most likely M parameter is 7.4+*q, indicating a large 
departure from the Pristine Model. The fraction of stars 
formed from gas that fell into the system is (M— 1) /M, or 



~ 87% in this case. The Extra Gas Model does very well 
at matching the symmetric peak, but the narrowness of 
the peak demands a large M, which in turn causes an un- 
derestimate of the frequency of low-metallicity stars and 
an overestimate of the frequency of high-metallicity stars. 
If we have overestimated the measurement uncertainties 
on [Fe/H], then the model's intrinsic, unconvolved peak 
could be wider, and M could be smaller. That solution 
would fit the metal-poor tail and the steep metal-rich 
slope better. 

The reasonable fit of the Extra Gas Model suggests 
the following extended star formation history for Fornax. 
Fornax began as a dark matter subhalo with a gas mass of 
~ 13% of its final stellar mass. Pristine, zero-metallicity 
gas fell in gradually. The first stars in Fornax formed 
from the initial gas. The most massive of these stars ex- 
ploded very quickly and enriched the interstellar medium 
(ISM) rapidly. The less massive, long-lived stars from 
early in Fornax's history still populate the metal-poor 
tail today. As more gas fell in, star formation continued 
from the enriched gas. Because the timescale for Type II 
SN metal enrichment is short, most of the stars in Fornax 
formed after several generations of massive stars already 
enriched the ISM. As a result, most Fornax stars oc- 
cupy the relatively metal-rich peak. This pi cture is qual- 
itative ly consistent with the conclusions of lOrban et al.l 
(2008). From color- magnitude diagrams, they deduced 
that 27% of Fornax stars are older than 10 Gyr, and 
33% are younger than 5 Gyr. The old stars enriched the 
ISM for the majority population — intermediate- age and 
younger. 

iBattaglia et al.l (|2006l ) also observed a spectroscopic 
(Ca triplet-based) MDF of Fornax. Their sam- 
ple extended beyond t he t idal radius at 1.2 deg 
(jlrwin fc Hatzidimitrioul 119951 ) . Our sample, which 
has a higher s patial density, extends to 16.3 arcmin. 
IBattaglia et all s MDF within 24 arcmin shows the same 
qualitative shape as our MDF. The peak of their MDF 
is ~ 0.15 dex lower than the peak we measure, but they 
showed that Fornax's [Fe/H] decreases beyond the ra- 
dial extent of our data. Therefore, we expect that their 
more extended sample would be more metal-poor. In 
fact, they argued for the importance of gas outflows and 
infall for the star formation history of Fornax based not 
on the MDF but on the radial metallicity gradient and 
non-equilibrium kinematics. 

3.2. Leo I 

Leo I is the second most luminous intact dSph that or- 
bits the MW. Accordingly, its peak [Fe/H] is lower than 
that of Fornax, with a median [Fe/H] of —1.42 and an 
intrinsic dispersion of .36 dex. In a Ca triplet study, 
IGullieuszik et al.l ()2009D found that the mean metal- 
licity is [M/H] = —1.2 with a dispersion of 0.2 dex. 
The shape of the MDF resembles that of Fornax. The 
MDF has a narrow peak and a long metal-poor tail. 
The most likely extra gas parameter is M = 7.2+J'g. 
The presumed increase in gas supply was about as in- 
tense for Leo I as for Fornax. Again, the Pristine and 
Pre-Enriched Models are not good descriptions of the 
MDF: lnL max (Extra Gas)/L max (Pristine) = 163.32 and 
In L max (Extra Gas) /L max (Pre-Enriched) = 81.87. In 
terms of the star formation history proposed for Fornax, 
it seems that Leo I encountered about the same history, 
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with a gradual increase in the gas supply that formed a 
small fraction of very metal-poor stars. 

Leo I has the largest distance (254 kpc.[Bcllazzin i~et al.l 
120041 ) among the "classical" dSphs. If ram pressure strip- 
ping of gas plays a r ole in the star formation history of 
dwar f galaxies (e.g., iLin fc Faberlll983l iMarcolini et al.1 
2003), then Leo I may have experienced less gas strip- 
ping thaji_Scuh>tor or Fornax. Its highl y elliptical or- 
bit (|Sohn et all l2007t iMateo et all [2008) allows Leo I 
to spend most of its time far from the MW center, al- 
lowing little time for strong tidal interactions. (When 
it does pass through the disk, however, it passes at 
high velocity.) In its lifetim e, Leo I experien ced a 
close periga lacticon just once (jMateo et ali [20081 or a 
few times ()Sohn et al.l I2007D . The possibly less in- 
tense interaction with the MW could have preserved 
enough gas in Leo I for it to appear more like Fornax 
th an Sculptor. The star formation histories deduced 
by lOrban et al.l (|2008l ) support that idea. Their mea- 
sured mea n ages are 6-8 Gyr fo r Leo I and Fornax stars. 
However, Smcckcr-Ha ne et al.l f|2009T ) interpreted their 
Hubble Space Telescope Advanced Camera for Surveys 
(HST/ACS) imaging as evidence that the dominant pop- 
ulation of Leo I is older than 10 Gyr. Our MDF gives 
only a qualitative suggestion that Leo I has experienced 
an extended star formation history, like Fornax, but the 
MDF alone is insufficient to quantify the mean stellar 
age. 

The close perigalactic passes probably caused Leo I to 
lose stars in addition to gas. The stars most susceptible 
to tidal stripping are the outermost ones, which tend to 
be older and more metal-poor (see Sec. [5]). Therefore, 
the observed MDF represents only the stars that Leo I 
still retains but not the complete star formation history 
o f the galaxy. 

iBosler et all (|2007| ) also measured the MDF for Leo I. 
They used their own calibration between the Ca triplet 
equivalent width and metallicity. They preferred to cal- 
ibrate to [Ca/H], but they also derived MDFs in terms 
of [Fe/H]. The average [Ca/H] of the ir MDF -1.3 4. The 
average [Fe/H] of our MDF is -1.43. IBosler et all 's sam- 
ple reached a radial extent of at least 20 arcmin. The 
maximum extent of our sample is 14.4 arcmin, and it be- 
comes sparse ly sampled beyond 8 arcmin. The shape of 
IBosler et alj 's MDF is slightly asymmetric, with a metal- 
poor tail. Our MDF is highly asymmetric. The sparser 
sampling of their MDF (102 stars) compared to ours (866 
stars) and the different measurement techniques may ex- 
plain the different appearances. 

3.3. Sculptor 

We ha ve alread y shown and analyzed the MDF for 
Sculptor ([Paper J ). We now place the MDF for Sculptor 
in the context of the MDFs for the other seven dSphs. 
Sculptor is the third most luminous intact dSph that 
orbits the MW. Its MDF might be expected to appear 
similar to Fornax and Leo I. However, Sculptor's MDF is 
unlike any of the other seven dwarfs. None of the three 
models adequately describe the MDF. 

First, the Pristine Model is too narrow to repro- 
duce the broad [Fe/H] distribution, including a possi- 
ble secondary peak at [Fe/H] = —2.1. Second, the 
Pre-Enriched Model might have been able to attain a 
better- fitting shape with a high enough [Fe/H] . How- 



ever, star S1020549 at [Fe/H] = -3.87 ± 0.21, con- 
firmed with high-resolution spectroscopy bv lFrebel et al.l 
(|2010a|) . demands a low initial enrichment. The up- 
per limit on the initial enrichment that we derive is 
[Fe/H] < -3.67. Third, the most likely Extra Gas 

Model, with M — 1.3^° i, is nearly identical to the Pris- 
tine model. Increasing M would only narrow the MDF 
further. 

A more proper nucleosynthetic treatment of iron would 
also better reproduce the MDF. Our analytic models 
assume instantaneous recycling, meaning that they do 
not incorporate a delay time between the births of stars 
and the return of enriched material into the ISM. In 
IPaper IVI (IKirbv et al.ll20loih . we show that a numerical 
model that relaxes the instantaneous recycling approxi- 
mation results in a model MDF that better fits Sculptor's 
observed MDF. However, even the more sophisticated 
model does not reproduce the apparent bimodality in the 
observed MDF (two peaks at [Fe/H] = -2.1 and -1.3). 

The kinematic distribution of Sculptor's stars may pro- 
vide s ome insight on the bimodal MDF. Tolsto y at al.l 
( 2004) found that a two-component model best describes 
Sculptor. The stars separate into a centrally concen- 
trated, metal-rich ([Fe/H] > —1.7), kinematically cold 
(<r„ = 7 km/s) component and an extended, metal-poor 
([Fe/H] < —1.7), kinematically warm (a v = 11 km/s) 
component. A two-component model would do a much 
better job of representing our observed MDF of Sculp- 
tor. Combining any two of the three chemical evolution 
models with different peaks in [Fe/H] would generate a 
broad [Fe/H] distribution. However, we reserve a two- 
component analysis for a study that includes kinematic 
data. 

The s hape of the Sculptor MDF is very different from 
Fornax. [Dekcl & Silk (1986) suggested that galactic out- 
flows play a large role in dSph formation. If Sculptor is 
less massive than Fornax, then winds from supernovae 
could have rapidly depleted Sculptor of the enriched gas 
necessary to create the more metal-rich stars. Therefore, 
the smooth shape of Sculptor's MDF compared to the 
peaked shape of Fornax's MDF may indicate that galac- 
tic outflows were more important than an increase in the 
gas supply in Sculptor's history, whereas the reverse was 
true for Fornax. 

3.4. Leo II 

Leo II, Sextans, and Carina are the next most lu- 
minous dSphs that orbit the MW. Their luminosities 
ar e nearly indistinguishabl e with in the error bars given 
by llrwin fc Hatzidimitrioul (|1995l ). Leo II continues the 
trend established by Fornax and Leo I. It has a slightly 
lower median [Fe/H] (—1.59) than that of Leo I. In accor- 
dance with its lower luminosity, its extra gas parameter 
M = 3.lt°;5 is also lower than that of Leo I. The Pris- 
tine and Pre-Enriched models are still not good fits, but 
they are better representations of the MDF than for For- 
nax and Leo I, where the Extra Gas Model departs more 
severely from the Pristine Model (M > 7). 

None of the models can reproduce the steep slope 
on the metal-rich side of the peak of Leo II. One pos- 
sible solution to the abrupt drop in the frequency of 
metal-rich stars is a terminal wind. For example, a 
Closed Box Model that is truncated at some metallic- 
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ity ([Fe/H] ~ —1.3 in thi s case) m i ght p roduce an MDF 
similar to that of Leo II. iWinnic k (2003), among others, 
considered such a model in her description of the MDFs 
o f Sculptor, Draco, a nd Ursa Minor. 

IQrban et al.l (|2008D reported similar mass- weighted av- 
erage stellar ages for Fornax (7.4 Gyr), Leo I (6.4 Gyr), 
and Leo II (8.8 Gyr). They listed the mean stellar ages 
for Sculptor, Sex tans, Ursa Minor, an d Draco as larger 
than 10 Gyr, and lMartin et al.l ([2008b) gave a mean age 
larger than 10 Gyr for Canes Venatici I. It is not surpris- 
ing, then, that Fornax, Leo I, and Leo II — the younger 
dSphs — have similar MDFs. Furthermore, their younger 
mean ages may have allowed them extra time to accrete 
gas while they were forming stars. The increased supply 
of gas would have caused their peaks to be narrower than 
can be explained by a Leaky Box Model and would ex- 
plain why the Extra Gas Model fits the best of the three 
chemical evolution models that we consider. 

Fin ally, the distance to Leo II is nearly as great as 
Leo I. ISieeel et al.l (1201(1 gave a distance of 219 kpc. Al- 
though the orbit of Leo II is unknown, it could spend 
much of its time far outside of the range of disruptive 
gravitational interaction with the MW. It could have 
been spared gas stripping, which other dSphs like Sculp- 
tor or Sextans may have encountered. However, if Leo II 
does spend most of its time in a low density region of 
the Local Group, then it would likely not encounter the 
gas reservoir required to explain our interpretation of the 
MDF as indicative o f an increase in the gas reservoir. 

iBosler et all (|2007t ) observed the MDF of Leo II in ad- 
dition to Leo I. They measured a mean [Ca/H] of —1.65, 
whereas we measure a mean [Fe/H] of —1.62. The radial 
extent of their survey was nearly the same as our survey. 
The shape of the Ca triplet MDF is at least qualitatively 
similar to our MDF. Both sho w an abrupt fa ll-off in the 
frequency of metal-rich stars. IBosler et al.l interpreted 
the absence of metal-rich stars as evidence for rapid gas 
loss. Indeed, the effective metallicity yield that we mea- 
sure (0.030 ± 0.002 Zq) is too low to be explained by 
completely retained SN ejecta. The galaxy must have 
lost some gas, but the shape of the MDF also mandates 
that it accreted low-metallicity gas. 

3.5. Sextans 

Despite its similar luminosity to Leo II, Sextans dis- 
plays a differently shaped MDF. It is more symmet- 
ric, with a shallower slope on the metal-rich side of the 
peak. Most strikingly, the median [Fe/H] for Sextans 
(—2.00) is significantly lower than the median [Fe/H] 
for Leo II (—1.59). As a result, the most likely yield 
for the Extra Gas Model is much lower in Sextans 
(0.015 ±0.001 Z Q ) than in Leo II (0.030 ±0.002 Z Q ) de- 
spite Sextans's similar extra gas parameter, M = 3.0j^'g. 
However, the Extra Gas Model may not be the best de- 
scription of the MDF. Quantitatively, neither the Pre- 
Enriched nor the Extra Gas Model is significantly pre- 
ferred: In L max (Pre-Enriched)/L max (Extra Gas) = 0.81. 

Sex tans has a very large tidal r adius (160 ± 50 ar- 
cmin. flrwin fc Hatzidi mitrioul 119951 ) . which means that 
a complete MDF req uires extensive radial sampling. 
iBattaglia et al.l (|2010T ) observed the Ca triplet-based 
MDF of Sextans to very large distance from the dSph 
center. Their sample extended to 1.8° whereas our sam- 



ple extends to only 21.4 arcmin. Their more complete 
MDF is more metal-poor than ours because they de- 
tected a radial gradient of —0.33 dex kpc -1 . Our sam- 
ple is too centr ally concentrated to detect a gradient. 
IBattaglia et al.l also found that the shape of the MDF 
changes as a function of radius. Sextans appears to 
have two metallicity populations: a metal-rich popula- 
tion within 0.8° and a metal-poor population beyond 
0.8°. Therefore, our results should be interpreted as ap- 
plicable only to the innermost population, which presum- 
ably formed more recently than the outer population. 

3.6. Draco 

Draco, Canes Venatici I, and Ursa Minor form the next 
group of dSphs with indistinguishable luminosities. The 
MDFs of these three dSphs and Sextans are broadly re- 
lated. Their mean [Fe/H] values all lie between —2.2 and 
— 1.9, and their MDFs seem to be more symmetric than 
the more luminous dSphs. 

The MDF of Draco is similar to that of Sextans. 
Most of the shape parameters (mean, median, disper- 
sion) for the two MDFs are nearly identical. In partic- 
ular, the observed MDF is more peaked than the Pris- 
tine Model, and it has a metal-poor tail. Formally, the 
Extra Gas Model is a better fit to Draco than the other 
models (In L max (Extra Gas) /L max (Pristine) = 20.34 and 
lnL max (Extra Gas)/L max (Pre-Enriched) = 2.84). On 
the other hand, the Pre-Enriched Model fits Sextans 
slightly better than the Extra Gas Model. However, 
Draco is better sampled than Sextans, with more than 
twice as many stars with measurements of [Fe/H]. As a 
result, the most likely chemical evolution parameters for 
D raco are more s ecure than for Sextans. 

IWinnick! (|2003l ) also observed and modeled the MDFs 
for Draco, Sculptor, and Ursa Minor. Her analytic mod- 
els included abrupt or continuous gas loss with initially 
pristine or pre-enriched gas. These single-component 
models did not fit any of the three MDFs very well, 
but a two-component model for Draco did work well. 
The two components were two Leaky Box Models with 
pre-enrichment. The models reproduce d the Dra co MDF 
shape very well. However, the shape of lWinnickl 's Draco 
MDF was different from ours. Her MDF had two peaks. 
The different shape results from different measurement 
techniques and diffe rent rad ial sampling (a maximum 
of ~ 30 arcmin for iWinnickl 's and 60 arcmin for our 
sample). Nonetheless, the two-component experiment 
demonstrates that an accurate description of some MDFs 
(especially Sculptor) may require multi-component mod- 
els. 

3.7. Canes Venatici I 

The MDF of Canes Venatici I resembles a normal dis- 
tribution more than any of the other seven dwarfs. Al- 
though it has a slight metal-poor tail, it is nearly sym- 
metric. Canes Venatici I shows the least preference 
for the Extra Gas Model of all the dSphs in Fig. [Tj 
except for Sculptor. The most likely extra gas pa- 
rameter is M — 1.6^Q 3, and the preference of the 
Extra Gas Model over the Pristine Model is slight 
(In I/ max (Extra Gas) /L max (Pristine) = 0.26). 

Overall, the MDF of Canes Venatici I fits a Leaky 
Box Model best of all the dSphs shown here. It is 
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the only dSph in our sample that fits the Pre-Enriched 
Model significantly better than the Extra Gas Model 
(In L max (Pre-Enriched) /L may ( Extra Gas) = 4.28 ). The 
dSph's large distance (210 kpc. lKuehn et al.ll2008l ) would 
make it less susceptible to gas accretion from the MW 
than the other dSphs in its luminosity class, Draco and 
Ursa Minor. Continuous gas outflow — from SN winds, 
for example — would not drastically affect the shape of 
the MDF. It would instead decrease the effective yield p 
and cause a low peak [Fe/H], which is observed. 

3.8. Ursa Minor 

The lumi nosity of Ursa Mino r is within a factor of two 
of Sextans. lOrban et al.l (|2008l ) derived identical star for- 
mation histories for both dSphs: no stars younger than 
12 Gyr. Ursa Minor, with the slightly lower luminos- 
ity, has a correspondingly lower median [Fe/H]: —2.13 
compared to —2.00. The most likely Extra Gas Model 
indicates intense gas inflow over the lifetime of star for- 
mation, despite the lack of an obvious G dwarf problem. 
The M parameter for the Extra Gas Model is 9.llg Q for 

Ursa Minor compared to 3.0]to'g for Sextans. 

None of the models accurately reproduce the sud- 
den absence of m etal-poor stars at [Fc/H] < —3. 
iHelmi et ail (|2006l) invoked the Pre-Enriched Model 
as a possible explanation for the apparent dearth of 
metal-poor stars in Fornax, Sculptor, Carina, and 
Sextans. Since then, extremel y metal-poo r stars 
have been discovered i n dSphs dKirbv et al.l 12008b, 
20091 iGeha et all [200l iCohen fc Huangl 1 20091 12010: 
Frebel et al.l l2010allbl; iSimon et al.ll2010allbl: iNorris et al 
2010allbl : iStarkenburg et al.1 120101 : iTafelmever et al 
20101 ) . These discoveries do not preclude pre-enrichment 



in all dSphs. In fact, the most likely initial metal- 
licity for the Pre-Enriched Model for Ursa Minor is 
[Fe/H ] n = -2.9ll^. However, ICohen fc Huangl 
(2010) have discovered two stars with [Fe/H] < —3 
in Ursa Minor using high-resolution spectra. Our 
sample includes seven such stars, including one with 
[Fe/H] = -3.62 ±0.35. 

On closer inspection, Ursa Minor is an outlier from its 
luminosity class. Whereas all of the other dSphs show 
at least a hint of a metal-poor tail, Ursa Minor shows 
a metal-rich tail. Ursa Minor's large negative radial ve- 
locity (—242.7 km s -1 ) rules out significant contamina- 
tion from metal-rich Galactic stars. The Pristine Model 
can explain neither the absence of a metal-poor tail nor 
the existence of a metal-rich tail. The most likely Pre- 
Enriched Model is too symmetric. The model with the 
sharpest peak is the Extra Gas Model with a very large 
M and a very low yield. The most likely parameters are 
M = 9.1±tt andp = 0.009±0.001 Z e , the most extreme 
of all eight dSphs. 

4. NUMERICAL CHEMICAL EVOLUTION MODELS 

We compare the observed MDFs of six galaxies (Leo I, 
Sculptor, Leo II, Sextans, Draco, and Ursa Minor) to pre- 
dictions of detailed numerical models of chemical evolu- 
tion in addition to the analytic models. Figure Q] includes 
the MDFs from the numerical models as dotted lines. 
The model MDFs have been convolved with the same 
function as the analytic models to approximate obser- 
vational uncertainty (Sec. [3]). In the adopted numerical 



models (|Lanfranchi fe Matteuccl2003ll2004l ). the evolu- 
tion of each galaxy is mainly controlled by the assump- 
tions regarding the star formation (SF) history, SF ef- 
ficiency, and galactic wind. All models adopt an infall 
of pristine gas, low star formation rate (SFR), and high 
galactic wind efficiency. The low SFR and the high effi- 
ciency of the wind give rise to a peak in the MDFs at low 
[Fc/H] (approximately between —1.6 to —2.0) whereas 
the long infall timescale of pristine gas allows the models 
to form a low number of metal-poor stars, similar to the 
observed frequency. 

In contrast to the analytic models, the numerical mod- 
els were not adjusted to fit the present data. The pre- 
dictions are the same as in previous papers, in which 
the models were adjusted to better match different ob- 
servational data. The numerical models f or Ursa Minor 
and Draco ([Lanfranchi fc Matteuccill2007l ) were adjusted 
to approximate the observed MDFs bas ed on photo- 
metric metallicities (lBella zz ini et al.ll2002f). whereas the 
models for Leo I and II ( Lanfranchi fc Matteuccil 120101) 
were calibrated to mat ch the [Fe/H] distri b ution inferred 
from Ca triplet lines dKoch et al.l l2007allbl : iBosler et all 
120071 iGullieuszik et al.ll2009f ). ~The predictions for Sex- 
tans and Sculptor are true predictions; they have not 
been calibrated to match any observational data be- 
cause they were published be fore data were available 
(|Lanfranchi fc Matteuccil 120041) . We can use the com- 
parisons to infer what modifications should be made in 
the models, especially those regarding SF and wind effi- 
ciencies. 

The differences in shape between the predicted and 
observed MDFs may be described in terms of how the 
modeled SF history could be changed to achieve a better 
fit. Higher SFRs seem to be necessary in the case of Leo I 
and Leo II because these two models exhibit MDF peaks 
~ 0.4 dex lower than observed. The model of Leo I is 
characterized by two long episodes of SF at 14 Gyr and 
9 Gyr ago, lasting 5 Gyr and 7 Gyr respectively, with a 
low efficiency (y ~ 0.6 Gyr -1 ) and by the occurrence of a 
very intense galactic wind with a rate 9 times higher than 
the SFR (wi — 9). On the other hand, the Leo II model 
adopts lower SF and wind efficiencies (v ~ 0.3 Gyr - and 
Wi = 8) and just one long episode of SF at 14 Gyr ago, 
lasting 7 Gyr. Modifying the duration and epoch of the 
SF episodes would not change the predictions substan- 
tially, whereas the wind efficiency influences the position 
of the peak in the MDFs and, most significantly, the rel- 
ative number of metal-rich stars. Because the shapes of 
the predicted MDFs are similar to the observations, SF 
efficiency is the only parameter that requires adjustment. 
Increasing v might lead to a better fit in both cases. 

For Ursa Minor and Draco, the predicted MDFs are 
0.3-0.5 dex more metal-rich than the observed MDFs. 
These two galaxies are characterized by short, older pe- 
riods of SF (4 Gyr and 3 Gyr, respectively) compared 
to the other dSphs and by the lowest SF efficiencies 
[y ~ 0.1 Gyr -1 and v ~ 0.05 Gyr -1 , respectively) among 
the six models analyzed. It seems that these values need 
to be further decreased to match the observations, espe- 
cially for Ursa Minor. This galaxy also exhibits a more 
extended high metallicity tail, which could be reproduced 
by a lower galactic wind efficiency. To prevent the subse- 
quent increase in the peak [Fe/H], the SFR should also be 
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decreased. In contrast to Ursa Minor, the prediction for 
Draco fits the number of metal-poor stars well but over- 
estimates the peak of the MDF and the number of high 
metallicity stars. The discrepancies are probably conse- 
quences of low wind efficiency. A higher wind efficiency 
would decrease the SFR after the onset of the wind, low- 
ering the number of metal-rich stars born. Such a model 
would better match the data by creating an MDF with 
a peak at lower [Fe/H] and with fewer metal-rich stars. 

The predicted MDF of Sextans reproduces the ob- 
served data well. The model adopts a long episode of SF 
(longer than 4 Gyr) with low rates and an intense galac- 
tic wind {v ~ 0.08 Gyr~ x and Wi = 9), giving rise to a 
main stellar population with low [Fe/H] (~ —1.8 dex). 
There seems to be an underprediction of the frequency 
of metal-poor stars, probably due to the extended infall 
of gas and to the prolonged SF. 

For Sculptor, as in the case of analytic models, the 
prediction does not fit the observed MDF. The ob- 
servations show a distribution much broader than the 
predictions. As discussed in Sec. 13. 3| the cause of 
the broad MDF could be two different stellar popu- 
lation with different m etallicities and spatial extents. 
lLanfranchi fc Matteuccf s model cannot separate the two 
different populations because it uses only one zone. In- 
stead, the different populations may be the result of dif- 
ferent SFRs in the central and in the outer regions, per- 
haps due to different gas densities. The result would be 
different chemical enrichment and different mean metal- 
licities in each region. We explore this idea more gener- 
ally in the next section. 

5. RADIAL METALLICITY DISTRIBUTIONS 

Galaxies often sh ow radial metallicity gradients (e.g., 
iMehlert et al.l l2~003j ). We discuss three processes that 
may be responsible for gradients. (1) The star f ormation 
rate i n a galaxy depe nds on the gas density (|Schmidtl 
Il959t rKennicuttiri983D . and the gas density increases to- 
ward the bottom of the gravitational potential well in 
the center of the galaxy. Therefore, the center of the 
galaxy may experience the highest star formation inten- 
sity and consequently may show the highest mean metal- 
licity. (2) MW satellite galaxi es may lose gas th rough 
tidal or ram pressure stripping (|Lin fe Faberl ll98~3t) from 
the MW. The gas that lies far from the center is more 
loosely bound to the dSph than gas that lies at the cen- 
ter. Therefore, gas leaves the dSph from the outside in. 
For dSphs affected by gas loss, later episodes of star for- 
mation occur closer to the center of the galaxy. Because 
later episodes of star formation occur from more metal- 
rich gas, gas loss creates a stellar metallicity gradient, 
with the center of the dSph being more metal-rich. This 
scenario requires star formation to occur during the pro- 
cess of gas loss. (3) The angular momentum of accreted, 
low-metallicity gas will not allow the gas to reach the 
center of the galaxy. If stars form from this gas, they 
will be metal-poor and mostly confined to large radius. 
Therefore, gas accretion can also generate radial metal- 
licity gradients. 

Figure [2] shows the radial distributions of [Fe/H]. 
The x-axis gives the distance from the dSph center 
in units of the core radius from a King profile fit 
(|Irwin fc Hatzidimitrioul I1995T) or half -light radius for 
Canes Venatici I ([Martin et al.l l2008bf ). The lines are 




-1 \* 

-2 

-3 
-4 




.3 



o "2 

"4 

-1 

-2 

-3 
-4 

-1 

-2 

1 

-1 

-2 

:1 

-1 

-2 

-3 
-4 



Fornax 

(-0.02 ± 0.02) dex per r c 
i i i 1 1 



*• 'Leo I 

(-0.11 ±0.01) dex per r c 
I I I ' ' 



***^r*V Sculptor 
(-0.18 ±0.02) dex per r c 



• Leo II 
(-0.21 ±0.01) dex per r c 
i i I'' 



Sextans 

(+0.06 ±0.05) dex per r c 
i i i''' 



Draco 

(-0.12 ±0.02) dex per r c 
i i i''' 



Canes Venatici I 
(-0.04 ± 0.03) dex per r h 
I I I''' 



Ursa Minor 
(-0.04 ±0.03) dex per r c 



12 3 4 

r/r c 

Fig. 2. — Radial distributions of [Fe/H] as a function of rad ius 
in units of the core radius (r c , Irwin & Hatzidimitriou 1995) or 
half-light radius (r^, Martin ct al. 2008b) in the case of Canes Ve- 
natici I. The dSphs are arranged in decreasing order of luminosity. 
The red line is the least squares linear fit, whose slope is given 
in each panel. The point sizes are inversely proportional to the 
uncertainty in [Fe/H]. 

the least squares linear fits. Table Q] gives the slopes of 
these lines in terms of angular distance from the dSph 
center in degrees, projected distance from the dSph in 
kpc, and projected distance in core radii. Four dSphs, 
Leo I, Sculptor, Leo II, and Draco, show significantly 
negative slopes. The slopes of the other four dSphs are 
consistent with zero or nearly consistent with zero. 

The eight dSphs in our sample are insulated against 
Galactic contamination by their Galactic latitudes (\b\ > 
40°), radial velocities (|i>heiio| > 200 km s _1 ), or both. 
Based on the velocity distribution of stars excluded by ra- 
dial velocity, we estimate that fewer than 5% of the stars 
are Galactic contaminants even in the worst cases (Sex- 
tans and Draco). However, so me Galactic ha l o star s may 
still contaminate our samples. iSchorck et al.l (|2009f) show 
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TABLE 1 

Metallicity Gradients 



dSph 


Distance 




r c 


d[Fe/U]/d9 


d[Fe/U]/dr 


d[Fc/R]/d(r/r c 




(kpc) 


(arcmin) 


(pc) 


(dex deg — 1 ) 


(dex kpc -1 ) 


(dex) 


Fornax 


139 


13.7 


550 


-0.07 ±0.09 


-0.03 ±0.04 


-0.02 ±0.02 


Leo I 


254 


3.3 


240 


-2.01 ± 0.10 


-0.45 ± 0.02 


-0.11 ±0.01 


Sculptor 


85 


5.8 


140 


-1.86 ±0.16 


-1.24 ±0.11 


-0.18 ±0.02 


Leo II 


219 


2.9 


180 


-4.26 ±0.31 


-1.11 ±0.08 


-0.21 ±0.01 


Sextans 


95 


16.6 


460 


±0.20 ±0.19 


+0.12 ±0.11 


±0.06 ± 0.05 


Draco 


69 


4.5 


90 


-0.73 ±0.13 


-0.45 ± 0.08 


-0.11 ±0.02 


Canes Venatici I 


210 


8.9 a 


540 a 


-0.48 ±0.36 


-0.13 ±0.10 


-0.07 ±0.05 


Ursa Minor 


93 


10.1 


270 


-0.21 ±0.19 


-0.18 ±0.16 


-0.06 ±0.05 



Refere nces. — Distance s adop ted from Rizz i et al. (2007) for Fornax, Bcllazzini et al.l 1 2004 ) 
for Leo I, IPietrzynski et al.l (p008h for Sculptor. I Sieeel et al. j2010t) for Leo II. ILee et al.l Miffl) 
for Sextans , Bcllazzini ct al. (2002) for Draco, IKuehn et alj 1 120081) for Canes Venatici I, and 
IMighell & Burke (1999) for Ursa Minor. Core radii adopted from Irwin & Hatzidimitriou (1995) 
except for Canes Venatici I, for which we adopt the half-light radius derived by IMartTneTlill 
(20083). 

a Half-light radius instead of core radius. 



that the local Galactic halo MDF is broad with a peak at 
[Fe/H] = —2.1 and a secondary peak at [Fe/H] = —1.1. 
We would not necessarily measure the correct metallicity 
for these stars because we assume the same distance mod- 
ulus for all stars in a given dSph sample. Halo stars could 
lie in front of or behind the dSphs. The distance modulus 
affects the metallicity measurements through the pho- 
tometric determination of surface gravity and partially 
photometric determination of effective temperature. We 
estimate that the intrinsic breadth of the halo MDF cou- 
pled with the increased measurement errors would cause 
halo stars to appear uniformly distributed in the range 
—3 < [Fe/H] < —1. Fornax and Leo I contain many 
stars at higher metallicities. Halo contamination would 
cause us to infer a more negative metallicity gradient for 
those two dSphs because the ratio of dSph to halo stars 
decreases with radius. However, the radial velocity distri- 
butions of stars in the Fornax and Leo I slitmasks suggest 
contamination at the level of 1% or less. Contamination 
in the other dSphs would tend to flatten the measured 
gradients and reduce their significance because the con- 
taminating population has roughly the same [Fe/H] dis- 
tribution as the average distribution for the entire dSph. 
However, we reiterate that the contamination is small. 

We might have expected that Fornax has a strong ra- 
dial gradient due to its longer star formation history 
(jOrban et al.1 12008). with successive generations becom- 
ing more and more centrally concentrated as the gas 
was depleted from the outskirts. Howe ver, the slope of 
the be st-fit line is consistent with zero. [Battaglia et al.1 
(2006) did find a radial metallicity gradient in Fornax. 
Their sample extended beyond the tidal radius whereas 
our sample reaches to only about the core radius. Within 
the radial bounds of our sample, Battag lia et all s mea- 
surements do not show evidence for a radial gradient, 
either. Therefore, our sample is too limited in angular 
extent to draw definitive conclusions on the presence of 
a radial metallicity gradient. 

Sextans, Canes Venatici I, and Ursa Minor also show 
insignificant slopes of [Fe/H] with radius within our cen- 
trally restricted samples. The lack of gradients possibly 
indicates short star formation durations. Gradients in 
[Fe/H] may occur because gas leaves the shallower po- 
tential in the dwarf galaxy's outskirts more easily than 



it leaves the center. Therefore, late star formation — 
from gas that had more time to be enriched — occurs 
only in the dSph's center. However, if the star forma- 
tion occurs over a period shorter than the gas redis- 
tribution timescale, then the dSph will show no [Fe/H] 
gradient. Our measurements of flat [Fe/H] radial distri- 
butions are consistent with a short star formatio n du- 
ration. Note, however, that iMarcolini et al.1 (|2008f ) pre- 
dicted that metallicity gradients are strongest when stars 
have been forming for 1 Gyr. After 1 Gyr, metals dis- 
persed by SN winds enrich the outer parts of the galaxy 
as much as the inner parts. Therefore, our interpretation 
of the shallow slopes is valid only if star formation in 
these dSphs spanned a period significantly shorter than 
1 Gyr. B oth photometry and detailed abundances (see 
IPaper IVt ) indicate that the SF durations were indeed 
t hat brief. 

iWinnick] (|2003l) has already arrived at some of these 
conclusions with her Ca triplet MDFs of Sculptor, Draco, 
and Ursa Minor. She detected radial metallicity gradi- 
ents in Sculptor and Draco, but not in Ursa Minor. She 
also found that the Leaky Box Model of chemical evolu- 
tion was a poor fit to the MDF of Ursa Minor, and she 
surmised that star formation in Ursa Minor was partic- 
ularly quick and efficient. 

6. METALLICITY TRENDS WITH LUMINOSITY 

So far, we have hinted at the possible role of total dwarf 
galaxy luminosity in determining its chemical evolution 
and metallicity distribution. In this section, we explicitly 
quantify trends of the MDF properties with luminosity. 
These properties include the mean metallicity, the in- 
trinsic width of the MDF, and the slopes of the radial 
gradients. 

Table [5] gives a summary of the MDFs for 15 dSphs: 
the 8 dSphs discuss ed in this ar t icle and the 7 additional 
dSphs discussed by Kir bv et al.l (|2008bl Leo T, Hercules, 
Ursa Major I and II, Leo IV, Canes Venatici II, and Coma 
Berenices). The MDFs of these 7 dSphs have been up- 
dated following the changes described in Papers I and 
II. The table shows the number of stars we analyze, lu- 
minosity, and different descriptions of the average [Fe/H] 
and width of the [Fe/H] distribution. The mean [Fe/H] is 
weighted by the inverse square of the errors. The MDF 
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TABLE 2 
Summary of dSph MDFs 



dSph 


iV a 


log(L/L©) 


([Fe/H]) b 




7° 


Median 


m.a.d. 


IQR e 


Skewness 


Kurtosis 




675 


7.3 ±0.1 


_n qq + n m 


0.36 


(0.31) 


— 1.01 


0.19 


0.37 


_i 4- n 09 


•3 _|_ n 1 q 


Leo I 


827 


6.7 ± 0.1 


— 1 43 + 01 


0.33 


(0.29) 


— 1.42 


0.18 


0.37 


— 1 47 + 09 


4 99 + 17 


Sr* ii 1 ntnr 


376 


6.4 ± 0.2 


_1 R» + 01 

l.UO _1_ W.L/.L 


0.48 


(0.46) 


— 1.67 


0.37 


0.75 


n 67 + 1 3 


25 + 25 


Leo II 


258 


5.8 ±0.1 


-1.62 ±0.01 


0.42 


(0.37) 


-1.59 


0.23 


0.51 


-1.11 ±0.15 


1.10 ±0.30 


Sextans 


141 


5.6 ±0.1 


-1.93 ±0.01 


0.48 


(0.39) 


-2.00 


0.29 


0.57 


-0.10 ±0.20 


0.47 ±0.41 


Draco 


298 


5.4 ±0.1 


-1.93 ±0.01 


0.47 


(0.36) 


-1.93 


0.26 


0.51 


-0.51 ±0.14 


0.73 ±0.28 


Canes Venatici I 


174 


5.4 ±0.1 


-1.98 ±0.01 


0.55 


(0.44) 


-1.98 


0.36 


0.71 


-0.26 ±0.18 


0.23 ±0.37 


Ursa Minor 


212 


5.3 ±0.1 


-2.13 ±0.01 


0.47 


(0.34) 


-2.13 


0.25 


0.50 


-0.03 ±0.17 


1.34 ±0.33 


Leo T 


18 


5.1 ±0.3 


-1.99 ±0.05 


0.52 


(0.43) 


-1.92 


0.31 


0.79 


-0.53 ±0.54 


-1.38 ±1.04 


Hercules 


21 


4.6 ±0.1 


-2.41 ±0.04 


0.64 


(0.56) 


-2.62 


0.46 


0.96 


0.51 ±0.50 


-0.94 ±0.97 


Ursa Major I 


28 


4.1 ±0.1 


-2.18 ±0.04 


0.64 


(0.59) 


-2.62 


0.30 


0.60 


0.61 ± 0.44 


-0.62 ±0.86 


Leo IV 


12 


3.9 ±0.2 


-2.54 ±0.07 


0.70 


(0.63) 


-2.35 


0.58 


0.86 


0.61 ±0.64 


-0.74 ± 1.23 


Canes Venatici II 


15 


3.9 ±0.2 


-2.21 ±0.05 


0.71 


(0.64) 


-2.68 


0.34 


0.62 


0.43 ± 0.58 


-0.37 ± 1.12 


Ursa Major II 


9 


3.6 ±0.2 


-2.47 ±0.06 


0.57 


(0.52) 


-2.41 


0.41 


0.92 


0.62 ± 0.72 


-0.69 ±1.40 


Coma Berenices 


18 


3.6 ±0.2 


-2.60 ±0.05 


0.40 


(0.26) 


-2.70 


0.29 


0.56 


-0.20 ±0.54 


-1.02 ±1.04 



Refe rences. — To derive the lumin osities of Fornax through Draco and Ursa Minor, we adopt the integrated U-band magni- 
tudes of Irwin & Hatzidimitriou (1995) and the distances given in Table [l] For the other galaxies, we adopt the luminosities of 
IMartin et al.l H2008al )! 

a Number of member stars, confirmed by radial velocity, with measured [Fe/H]. 

b Mean [Fe/H] weighted by the inverse square of estimated measurement uncertainties. 

c The number in parentheses is the intrinsic [Fc/H] spread corrected for measurement uncertainties. 

d Median absolute deviation. 

c Interquartile range. 

f Actually the excess kurtosis, or 3 less than the raw kurtosis. This quantifies the degree to which the distribution is more sharply 
peaked than a Gaussian. 



width er is given both in terms of the standard devia- 
tion uncorrected for measurement error as well as the 
width reduced by the amount of inflation caused by the 
estimated measurement uncertainties. We estimate the 
intrinsic spread cr( [Fe/H]) by solving the following equa- 
tion: 



1 N 

-y 



([Fe/H], t -([Fe/H])) 2 
A^(5[Fe/HU2+a([Fc/H]); 



= 1 



(8) 



The value of cr( [Fe/H]) for each dwarf is given in paren- 
theses in the column labeled a in Table [5] The last five 
columns show different shape parameters. The median 
identifies the peak of the MDF better than the mean. 
The median absolute deviation (m.a.d.) and interquar- 
tile range (IQR) are different measures of the width of the 
MDF. Skewness quantifies the asymmetry of the MDF, 
with negative values indicating a metal-poor tail. Kurto- 
sis quantifies by how much the MDF is peaked. Positive 
kurtosis indicates that the distribution is more sharply 
peaked than a Gaussian. 

6.1. Luminosity-Metallicity Relation 

The average metallicity of more luminous dwarf galax- 
ies is larger than for fainter dwarf galaxies. The relation 
between gas phase oxygen abundances a nd galaxy lumi- 
nosit y is particularly w e ll-stu died (e.g., iSkillman et al.l 
Il989t iVaduvescu et all 120071) . Others have stud- 
ied the more basic relation — the fundamental line — 
for d w arf galaxies (TPrada fc Burkertl I2002t IWoo et all 
120081) . IKirbv et al l (2008b) determined the luminosity- 
metallicity relation (LZR) based on medium resolution 
spectral synthesis of stars in eight faint dSphs combined 
with Ca triplet metallicity measurements for more lu- 
minous dSphs. Since then, our technique for measuring 
metallicities has been revised (Papers I and II), and we 



have calculated synthesis-based metallicities for the more 
luminous dSphs. The top panel of Fig. [3] shows the LZR 
for dwarf galaxies with metallicity measurements based 
only on spectral synthesis. 

The following equation describes the orthogonal re- 
gression fit 1 _a£COuirthrg_jOTCTrorsin both luminosity and 
[Fe/H] (JAkritas fc Bershadvlll996l ). where the errors are 
the standard deviations of the slope and intercept: 



([Fe/H]) = (-2.04 ±0.04) + (0.31 ±0.04) log 



( L t 



V10 5 L G 



(9) 

The linear Pearson correlation coefficient for the data is 
0.95, indicating a highly significant correlation. 

By including [Fe/H] measurements for other dwarfs, 
we may refine the fit at the cost of losing the homo- 
geneity of the abundance analys is. We add Ca tr iplet- 
based m easurements for Ca rina (|Helmi et all 120061 ) and 
Bootes I (|Martin et al.l l2~007). The bottom panel of Fig. [3] 
shows the result. This LZR now includes all MW dwarfs 
less luminous than Sagit tarius except the le ast lumi- 
nous objects (Willman 1, iWillman et "all 120051: Segue 1 
iBelokurov eTall 120071: Segue 2, iB elokurov et all 120091: 
Boot es IL iWalsh et all 120071: Leo V, iBelokurov et ail 
120081 : Pi sces I. iWatkins et al.l 120091: and Pisces II and 
Segue 3, IBelokurov et al.l 120101 ) because they have only 
a few RGB stars, and their average metallicities are not 
well-determined. For the most part, even t heir luminosi- 
ties ar e uncertain by factors of two or more (|Martin et all 
l2008al K As expected, the addition of two galaxies to the 
existing 15 hardly changes the LZR: 



([Fc/H]) = (-2.02 ± 0.04) + (0.31 ± 0.04) log 



( L t 



\ 10 5 L, 



(10) 
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Fig. 3.— Top: The mean [Fe/H] of MW dSphs as a function 
of total luminosity. The dashed line is the weighted, orthogonal 
regression linear fit in log( L)-[Fe/H] space, accounting for the er- 
rors in both L and [Fc/H] ( Akritas & Bcrshadv 1996). The dot- 
ted lines are the rms dispersion of the residuals. The filled di- 
amonds represent measurements from this series of papers. The 
open di amonds represent the updated measurements of Kirby et al. 
(2008b), which were performed identically to those measurements 
presented here. Bottom: Same as the top panel, with two more 
galaxies that were not measured in the same way. Metallicity 
measurements for Bootes I and Carina are based on the equiva- 
le nt width of the Ca triplet. The dot-dashed line is the relation 
of Woo, Courteau, & Dekel (2008) from Local Group galaxies, in- 
cluding galaxies much more luminous than Fornax. 



The linear Pearson correlation coefficient for the data is 
0.93. 

A straight line may be an overly simplistic model to 
the luminosity-metallicity relation. The dwarfs with 
log(L/L ) > 5 seem to lie along a steeper line than the 
less luminous dwarfs. In order to better show that differ- 
e nce, the dot-dashed line in Fig. El is the best- fit relation 
of lWoo. Courteau. fc Dekel (|2008j) . They studied scaling 
relations among 41 luminous Local Group dwarf galaxies. 
They found the following relation between metal-fraction 
Z and stellar mass M* for dwarf ellipticals: 



log Z = -0.11 + (0.40 ± 0.05) log (j^) ■ (11) 

Assuming , as IWoo et all did, that [Fe/H] = log ( Z /Z Q ) 
and Z Q = 0.019 (jAnders fc Grevessel [1989') and that 
M*/L = 1.36 Mq/Lq (the average value for their 18 
dEs) , we replace Eq. [TT] with 



Fig. 4. — Top: The intrinsic spread in [Fe/H], accounting for 
measurement error, as a function of to tal dSph luminosity. The 
dashed line is the least-squares fit (Eq. 1130 . and the dotted lines 
represent the rms about the dashed line. Bottom: The logarithm 
of the linear metallicity spread as a functi on o f dSph luminosity. 
The dashed line is the least-squares fit (Eq. 1 141 ) . 



[Fe/H] = -2.06 + (0.40 ± 0.05) log (jjj^j ■ (12) 

This is the dot-dashed line in Fig. [3] It is an excellent 
fit to the luminous half of those 17 dwarfs. This is not 
su rprising bec ause many of those dwarfs were included 
in IWoo et all 's sample. However, the fit is not good to 
the dwarfs with log(L/L ) < 5. 

Despite the possible deviation from the LZR at low lu- 
minosities and low metallicities, the LZR is continuous — 
if not linear — from ultra-faint dSphs to massive ellipti- 
cal galaxies. iTremonti et al.l (|2004f ) demonstrated the 
tight correlation between gas phase metallicity and stel- 
lar mass over a wide range of masses. They deduced 
that low m ass galaxies pr eferentially lose metals to galac- 
tic winds (Larson 1974b). Their conclusion is also con- 
sistent with the stellar metallicity-stellar mass relation, 
such as derived fr om spectrophotometric indices (e.g., 
I Mendel et al.l [2009) . The continuity of the relation from 
massive ellipticals to galaxies with the luminosity of 
Coma Berenices (4000 L@) may indicate that the main 
variable that dictates the amount of metals galaxies lose 
is galaxy mass. 

6.2. Intrinsic [Fe/H] Spreads 



Metallicity Distributions in MW Dwarf Satellites 



13 



We expect the mean metallicity of a dSph to vary with 
luminosity because the presence of many stars implies a 
history of many SNe to enrich the gas. The more complex 
enrichment histories of the more luminous dSphs moti- 
vate us to examine how the width of the MDF varies with 
luminosity. 

Figure H] shows the trend of <r([Fe/H]) with dwarf 
galaxy luminosity. Intrinsic [Fe/H] spreads are larger in 
less luminous dwarfs. The least-squares fit is 



<r([Fe/H]) = (0.45 ± 0.03) - (0.06 ± 0.02) log I - L ' 



V10 5 L G 

(13) 

The linear Pearson correlation coefficient for the data is 
—0.57. The coefficient is nega tive because a( [Fe/H ]) is 
anticorrelated with luminosity. iNorris et al.l (|2010aD also 
quantified the luminosity-metallicity spread-luminosity 
relation. They suggested that an increase in cr([Fe/H]) 
at low luminosity indicates inhomogeneous and stochas- 
tic chemical enrichment in the lowest mass galaxies. 

However, [Fe/H] is a logarithmic quantity. In order 
to better visualize the physical metallicity spreads, we 
recast the top panel of Fig. [4] in terms of the linear metal 
fraction: Z/Z Q = 10[ Fc / H l. The bottom panel of Figure H 
shows the logarithm of ct(Z/Zq) versus metallicity. The 
least-squares linear fit is 



loga(Z/Z Q ) = (-1.90±0.05)+(0.28±0.05)log 



( L t 



(14) 

The linear Pearson correlation coefficient for the data is 
0.85, which indicates that the luminosity- Z spread rela- 
tion is much more significant than the luminosity- [Fe/H] 
spread relation. 

Note that \oga(Z / Zq) is not the same as 
<r[\og(Z/ZQ)] « er([Fe/H]). It is possible to approxi- 
mate cr([Fe/H]) as a(Z)/Z. Instead, we translated all 
values of logarithmic [Fe/H] into linear Z = Z Q 10[ Fe / H l 
and recalculated (j(Z/Zq) in analogy to Eq. [5] We 
estimated SZ through standard error propagation: 
SZ = (Zlnl0)<5[Fe/H]. 

The linear representation of metallicity spread shows 
that the physical (linear) ranges of metallicity are some- 
what similar in many dwarfs. Nine of the 15 dwarf galax- 
ies cluster around log <j(Z/Zq) — —2. More interestingly, 
the four most luminous dwarfs, Fornax, Leo I, Sculptor, 
and Leo II, have \oga(Z/Z Q ) > —1.6. 

These trends might be explained by the differences in 
star formation duration across dwarfs. In chemical evo- 
lution models, the width of the MDF depends, at least in 
part, on the effective yield p. As an example, the width 
of the MDF in the Pristine Model (Eq.© is a(Z) = V&P- 
For the four most luminous dSphs, the best-fit yields are 
p > 0.02 for all three chemical evolution models. For the 
next four dSphs in order of luminosity, the best-fit yields 
are p < 0.02 for all three models. The effective yield p 
encompasses gas outflow as well as SN iron yields. (See 
Table [3] for a list of the best-fit yields.) Therefore, the 
luminosity-metallicity spread relation may indicate that 
the more luminous dSphs more effectively retained their 
gas than lthe less luminous dwarf galaxies. As a result, 
the more luminous dSphs could have maintained star 



formation for a longer duration. The color-magnitude 
diagrams of Fornax, Leo I, and Leo II show that they 
did indeed e xperience star format i on more recently than 
10 Gyr (e.g.JMighefl fe Rich|[T99l l Buonanno et al| [T99l 
ISmecker-Hane et al.1120091 ). unlike th e other dS phs. Our 
results from the [a/Fe] distributions ([Paper PvT ) also sup- 
port star formation durations longer than 1 Gyr for the 
four more luminous dSphs and shorter than 1 Gyr for 
the four less luminous dSphs. 

6.3. Radial Gradients 

The radial gradients of the eight dSphs do not sepa- 
rate neatly into more and less luminous categories. We 
might reasonably guess that the more luminous galaxies, 
which experienced more prolonged star formation, would 
show steeper radial gradients. However, the most lumi- 
nous dSph in our sample, Fornax, does not show a radial 
gradient within the bounds of our data. Furthermore, 
Sculptor and Leo II show very strong radial gradients 
even though they are near the middle of the luminos- 
ity range of our sample of dSphs. The ex planation may 
relate to Sculptor's kinematic complexity (iTolstoy et al.l 
[200l iBattaglia et al.l l2008at I Walker et all 120071 12009). 
Future studies may also reveal multiple kinematic popu- 
la tions in Leo II to acc ompany its strong radial gradient. 

Spolaor et al.l (|2009f) measured the radial gradients of 
more luminous (Mb < —16.8) galaxies in the Fornax and 
Virgo clusters. The most luminous galax y in our sam- 
ple is the Fornax dSph (Mb — —12.6). ISpolaor et al.l 
found that the magnitude of radial gradients decreases 
with decreasing luminosity until Mb ~ — 17, where the 
radial gradient vanishes. The trend with velocity dis- 
persion is stronger than with luminosity The gradients 
disappear by ~ 45 km s _1 , much larger than the veloc- 
ity dispersion of any galaxy in our sample. Therefore, it 
is notable that half of the galaxies in our sample (Leo I 
and II, Sculptor, and Draco) display radial gradients. We 
conclude that the radial gradients of MW dSphs do not 
obey the tight relation seen for more luminous galaxies. 
Instead, we speculate that the particular star formation 
histories influenced by particular interactions with the 
MW determine the presence or absence of radial metal- 
licity gradients. 

At odds with ISpolaor et ail . Koleva et all (|2009aD 
found no trend between radial metallicity gradients and 
galaxy mass for galax i es in th e Fornax cluster and nearby 
groups. IKoleva et al.l (|2009b[ ) suggested possible reasons 
for the discrepancy, but they did not find any explana- 
tions satisfactory If IKoleva et all s result holds, then it 
would not be surprising that our sample also does not 
show a trend between gradients and luminosity. Regard- 
less, our samples are too centrally concentrated to probe 
the full extent of the radia l gradients for most of the 
dSphs. Other sam ples fe.g.. IBattaglia et al.ll2006l 120101: 
Wal ker et al.ll2009T ) have better addressed this issue. 

7. CHEMICAL EVOLUTION TRENDS WITH GALAXY 
PROPERTIES 

The chemical evolution model fits more directly relate 
to the star formation histories of the dSphs than the 
bulk metallicity properties. However, the derived quan- 
tities (p for the Pristine Model, p and [Fe/H] for the 
Pre-Enriched Model, and p and M for the Extra Gas 
Model) are not direct observables. In calculating them, 
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TABLE 3 

Chemical Evolution Models 





Pristine Model 


Pre-Enriched Model 


Extra Gas Model 


i L m ,v ( Pre-Enriched) 

|jl I1JctA ' L 

L max (Pristine) 


i LmnvfExtra Gas) 
L n „ (Pristine) 


dSph 


P (Ze) 


p(z®) 


[Fe/H] 


P (Zq) 


M 


Fornax 


0.118 ±0.005 


0.096 ± 0.005 


- 2.12 ±0.07 


0.122 ±0.004 


74 +i.2 


49.98 


125.62 


Leo I 


0.044 ± 0.002 


0.033 ± 0.002 


9 oc+0.05 
- a -°°-0.06 


0.045 ± 0.001 


7 n + 1.0 
'• z -0.9 


81.45 


163.32 


Sculptor 


0.029 ±0.002 


0.028 ±0.002 


< -3.67+ 33 


0.029 ± 0.002 


i q+0.2 
L - 6 -0.1 


-0.11 


1.59 


Leo II 


0.029 ± 0.002 


0.025 ± 0.002 


- 2 94+ 011 


0.030 ± 0.002 


o 1+0.6 

^•-'■-o.s 


12.67 


21.08 


Sextans 


0.016 ±0.002 


0.014 ±0.002 


-3.20^'il 


0.015 ±0.001 


3 0+ 1 ' 2 


6.43 


5.62 


Draco 


0.016 ±0.001 


0.013 ±0.001 


-3.05+^? 


0.015 ±0.001 


4 0+ 1 - 1 
4 - u -o.s 


17.51 


20.34 


Can. Ven. I 


0.018 ±0.002 


0.016 ±0.002 


o 7 +0.20 

001 -0.37 


0.017 ±0.002 


1 6+°' 5 


4.55 


0.26 


Ursa Minor 


0.011 ±0.001 


0.007 ±0.001 


_ 2 ni+0.09 
z ' yi -0.10 


0.009 ± 0.001 


a 1+4-4 
a - 1 -3.0 


22.10 


21.63 




Fig. 5. — Parameters of the best-fitting chemical evolution models (Table [3j as a function of galaxy properties. The model parameters 
are the effective metal yield (p) in units of the solar metal fraction, the initial metall icity ([Fe/H] n ) in the Pr e-Enriched Model, and the 
infall parameter (M) in the Extra Gas Model. The galaxy properties arc luminosity (Irwin & Hatzidimitriou 1995; Martin ct al. 2008b), 
line-of-sight velocity dispersion, two-dimensional half-light radius (both from Wolf ct al. 2010, and references therein), and Galactocentric 
distance. The error bars represent the asymmetric 68.3% confidence interval from Monte Carlo trials. The [Fe/H] value for Sculptor is an 
upper limit. 



we have assumed that the models are good descriptions 
of the MDFs. Although we have estimated the relative 
goodness of fit between the Simple and Extra Gas Mod- 
els, we have not estimated the absolute goodness of fit. 
Thus, the results of this section should be viewed as more 
directly relevant to the star formation histories whereas 
the results of the previous section should be viewed as 



more directly observable and therefore more confident. 

Table [3] presents the most likely chemical evolution 
model parameters for each dSph along with the ratios of 
the maximum likelihoods. The two-sided uncertainties 
represent the 68.3% confidence interval (see Sec. [3]). Sin- 
gle error bars are given for the yields because the upper 
error bars are equal to the lower error bars in all cases. 
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Many of the Monte Carlo trials for the Pre-Enriched 
Model of Sculptor reached to very low values of [Fe/H] Q . 
The formal value for the lower error bar is 14 dex. There- 
fore, we treat the calculated value of [Fe/H] as an upper 
limit. 

Fig. [5] shows the most likely parameters from Table |3] 
plotted against the luminosity, line-of-sight velocity dis- 
persion, projected half-light radius, and Galactocentric 
distance. The most obvious observation is that the yield 
p monotonically increases with dSph luminosity for the 
six most luminous dSphs. That increase is a reflection of 
the increasing average metallicity and increasing metal- 
licity spread. For the remaining model parameters, the 
dSphs may be separated into two broad categories: more 
luminous, infall-dominated (Fornax, Leo I, and Leo II) 
and less luminous, outflow-dominated (Sextans, Ursa Mi- 
nor, Draco, and Canes Venatici I). 

The more luminous, infall-dominated dSphs are more 
consistent with the Extra Gas Model than the Pristine 
Model or the Pre-Enriched Model. Even so, the most 
likely initial metallicities for the Pre-Enriched Model are 
unreasonably large. It would be strange for Fornax and 
Leo I to have been born with a much more metal-rich 
ISM than other dSphs. There is no evidence that they 
formed their first stars long enough after the less lumi- 
nous dSphs for the intergalactic medium to have become 
so enriched. The most likely Extra Gas Models have 
M > 3, indicating that gas infall significantly affected 
star formation over the lifetime of the galaxy. 

The less luminous, outflow-dominated dSphs show sim- 
ilar, low effective yields (0.007 Z Q < p < 0.018 Z Q ) 
compared to the more luminous dSphs, regardless of the 
chemical evolution model considered. It is possible that 
the average SN yields of the least luminous dwarf galax- 
ies are anomalous because the IMF was stochasticall y 
sampled in tiny stellar systems (e.g.. iKoch et al.l f2008). 
A more likely explanation for the low values of p is that 
gas outflow reduced the effective yield below the value 
achieved by SN ejecta. The low masses of the less lu- 
minous dSphs rendered them unable to retain their gas. 
Gas flowed out of the galaxy from internal mechanisms, 
such as SN winds, and external mechanisms, such as ram 
pressure stripping. The outflows prevented the MDFs 
from achieving a high ([Fc/H]) and caused the MDFs to 
be more symmetric than the more luminous dSphs. 

However, a curious dynamical property of dSphs may 
undermine this star formation hypothesis. The dwarf 
galaxies of the Local Group less luminous than For- 
nax seem to inhabit dark matter halos of similar mass 
(lMatedll998t lPeharrubia et al.ll200llStrigari et~aLll2008t 
I Walker et all 120091 ) . Although their stellar masses span 
nearly five orders of magnitude, their dominant dark 
matter masses are the same. As a result, gas out- 
flow from the center of, for example, Canes Venatici I 
should be no stronger than from Fornax. However, the 
masses are poorly constrained beyond the half-light ra- 
dius. It is possible that dSph mass profiles diverge at 
large radii, allowing blown-out gas to return to the more 
luminous galaxies with possibly larger total mass. How- 
ever, dynamical tracers beyond 300 pc are sparse, es- 
pecially for the less luminous dSphs. Consequently, in- 
no vations in measuri ng the total mass of dSphs (such 
as iWolf et all's |2010 | anisotr o py-ind ependent mass esti- 
mates or lAmorisco k, EvansT s 120101 phase space models) 



will be necessary to better constrain the return fraction 
of blown-out gas. 

We have refrained from classifying Sculptor into one 
of these two broad categories because none of the star 
formation models is a good fit to its MDF. It is possible 
that the hierarchical assembly of two dSphs has brought 
together a su p erpos ition of distinct stellar populations. 
iTolstov et al.l (|2004[) recognized this possibility in their 
finding of two distinct kinematic and metallicity com- 
ponents in Sculptor. Hierarchical assembly almost cer- 
tainl y plays some role in the early form ation of dSphs 
(e.g.. IDiemand. Kuhlen. fc Madaull2007fl . However, this 
process does not dominate the shapes of the MDFs for 
all dSphs. For example, Fornax and Leo I do not show 
evidence of distinct metallicity peaks. If Fornax or Leo I 
accreted so many "sub-dwarfs" that the metallicity peaks 
are no longer distinct, then we would not expect the av- 
erage metallicities of Fornax and Leo I to be so high and 
the MDF widths to be as narrow as they are. Instead, 
we conclude that Sculptor is unique among these eight 
dSphs in showing the most obvious sign of the superpo- 
sition of two separate stellar populations. 

With one exception, the only parameter to show a 
trend with any galaxy property is the yield, which is cor- 
related with luminosity. The other parameters, [Fe/H] 
and M, do not show obvious trends with luminosity, 
and no parameter shows a trend with velocity dispersion, 
half-light radius, or Galactocentric distance. The excep- 
tion is that the dSphs with higher cr\ oa seem to require 
larger values of M, or more extra gas. It would seem that 
the greater gravitational potentials of the galaxies with 
higher o~i os attracted more external gas to power star for- 
mation. However, we caution against this interpretation 
because the mass profiles of dSphs are complex, and o~i os 
does not completely represent the ability of a dSph to 
attract additional gas. Instead, we regard the trend of 
M with ctios as tenuous at best. As an example, <7i os for 
Leo I and Sculptor are identical, yet they have highly 
discrepant values of M. 

The chemical evolution models we have considered are 
overly simplistic. In reality, the dSphs probably have 
complex star formation histories. A steady star forma- 
tion rate with instantaneous mixing and instantaneous 
recycling does not completely describe the MDF of any 
dSph. Even our Extra Gas Model assumes a contrived 
functional form of the gas increase (Eq. [3]). Some pro- 
cess, such as interaction with the MW, must cause gas to 
fall into the galaxy and trigger the star formation rate to 
increase one or many times over the dSph's lifetime. We 
anticipate that more realistic semi-analytic and hydrody- 
namical models will provide much better comparisons to 
the observed MDFs than the very simple analytic models 
we have considered. 

8. CONCLUSIONS 

We have analyzed the metallicity distributions for eight 
dwarf satellite galaxies of the Milky Way. We fit analytic 
chemical evolution models to the MDF of each galaxy. A 
Leaky Box Model starting from zero-metallicity gas does 
not faithfully describe any of the galaxies because it en- 
counters the same "G dwarf problem" that once com- 
plicated the inte rpretation of the Milky Way's met a llic- 
ity distribution (|van den Berghlll962t ISchmidtlll963D . A 
model with a fairly arbitrary prescription for an increase 
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in gas supply better describes the shape of the MDFs by 
allowing for a narrower peak and a longer metal-poor tail 
than the Leaky Box Model. Permitting a non-zero ini- 
tial metallicity (pre-enrichment) allows the shape of the 
Leaky Box Model to better fit the observed MDFs, but 
in no case except Canes Venatici I does the Pre-Enriched 
Model fit obviously better than the Extra Gas Model. In 
several cases, the Extra Gas Model fits much better than 
the Pre-Enriched Model. 

The shapes of the MDFs follow several trends with 
luminosity. The strongest trend is the luminosity- 
metallicity relation. Final dwarf galaxy luminosity can 
predict a dwarf galaxy's mean [Fe/H] to within 0.16 dex 
(lcr confidence interval). The luminosity also deter- 
mines the width of the MDF. However, the luminosity- 
metallicity spread relation is not as smooth as the 
luminosity-metallicity relation. Instead, luminosity sep- 
arates the metallicity spreads into high or low, with 
the four most luminous dSphs having large spreads and 
dSphs with the luminosity of Sextans or smaller having 
small spreads. 

We surmise that dSph luminosity is a good indicator 
of the ability to retain and accrete gas, despite t he find- 
ing o f a common central den sity for all dSphs (Mated 
[19981 iPenarrubia et~aT1 [20081: iStrigari et all 12008ft . For 
the more luminous dSphs, an increase in the gas reservoir 
during the star formation lifetime shapes the MDF and 
keeps the effective yield and therefore the mean metal- 
licity and metallicity spread high. The less luminous 
dSphs are less able to retain gas that leaves via super- 
nova winds or interaction with the Milky Way. Finally, 
all of the chemical evolution models we consider are much 
too narrow to explain the MDF for Sculptor. The previ- 
ous evidence for multipl e kinematic populations in Sculp- 
tor ([Tolstov et alJl2004t iBattaglia et al.ll2008aft suggests 
that Sculptor experienced at least two distinct, major 
episodes of star formation. We have not considered kine- 
matics in our analysis, but we recognize the value radial 
velocities add to abundance data. We plan to explore 
the relationship between velocity dispersion, [Fe/H], and 
[a/Fe] in the future. 

Stellar mass, not lumin osity, is lik e ly the independent 
variable, a concept that IWoo et all ((2008[ ) explored in 
detail. The conversion from luminosity to stellar mass 
involves the ages and metallicities of the component pop- 
ulations. Rather than complicate the observational data, 
we have chosen to present the dSph metallicity relations 
against luminosity because it is a direct observable. Most 
of the dSphs in our sample are ancient and metal-poor. 
As a result, the relation between their luminosities and 
stellar masses is one-to-one. Three dSphs — Fornax and 
Leo I and II — have younger populations, which result 
in higher luminosities at a given stellar mass. Nonethe- 
less, the relation between luminosity and stellar mass is 
roughly monotonic for the dSphs i n our sampl e , acco rd- 
ing to the stellar masses derived bv lWoo et al.l ()2008[ ) or 



lOrban et all (|2008h . Therefore, we conclude that lumi- 
nosity is a good proxy for the more fundamental param- 
eter, stellar mass, for these eight dSphs. 

The radial gradients of [Fe/H] do not obey a relation 
with total dSph luminosity. The slopes of the gradients 
for four dSphs are consistent with zero, and four more 
are significantly negative. The slopes for Sculptor and 
Leo II are at least as steep as —0.18 dex per core ra- 
dius. Negative slopes are to be expected for most galax- 
ies. Consequently, we have no satisfying explanation for 
the lack of a pattern for which dSphs happen to show 
radial gradients, even though we have characterized the 
star formation of the four less luminous dSphs as dom- 
inated by gas outflow. Samples more radially extended 
than ours might show metallicity gradients in all dSphs. 

Although we have assigned quantitative parameters for 
analytic chemical evolution models to each galaxy, we do 
not believe that any of the chemical evolution models 
are excellent fits to any dSph. In particular, the assump- 
tions of instantaneous mixing and recycling are inappro- 
priate for small galaxies with star formation lifetimes 
longer than the timescale for Type la supernova explo- 
sions. We expect that models that incorporate inhomo- 
geneous pockets of star formation (such as the models of 
iMarcolini et"al|[2006l [20081: iRevaz et al.ll2009ft and time- 
delayed ir on enhancement from Type la SNe (such as the 
models oflLanfranchi fc Matteuc ci 2004; Marc olini et al.l 
120061 120081 : IRevaz et al.ll2009ft will much better describe 
the dSph MDFs. Interested modelers wishing to com- 
pare their predictions to our observed MDFs may find t he 
complete catalog of [Fe/H] measurements in iPaperTH . 
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